Charge and spin dynamics in interacting quantum dots 



Janine Splettstoesser 1,2 , Michele Governale 3 , Jiirgen Konig 4 , and Markus Blittiker 5 
1 Institut fur Theoretische Physik A, RWTH Aachen University, D-52056 Aachen, Germany 
2 JARA - Fundamentals of Future Information Technology, Germany 
3 School of Physical and Chemical Sciences and MacDiarmid Institute for Advanced Materials and Nanotechnology, 
Victoria University of Wellington, PO Box 600, Wellington 6140, New Zealand 
4 Theoretische Physik, Universitdt Duisburg-Essen and CeNIDE, D-47048 Duisburg, Germany 
5 Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve 4, Switzerland 

(Dated: May 3, 2010) 

The transient response of a quantum dot with strong Coulomb interaction to a fast change in 
the gate potential, as well as the stationary ac-response to a slow harmonic variation in the gate 
potential are computed by means of a real-time diagrammatic expansion in the tunnel-coupling 
strength. We find that after a fast switching, the exponential relaxation behavior of charge and spin 
are governed by a single time constant each, which differ from each other due to Coulomb repulsion. 
We compare the response to a step potential with the RC time extracted from the ac response. 



PACS numbers: 73.23.-b,73.23.Hk 
I. INTRODUCTION 

Single electron charging of nanoscale devices driven by 
a time-dependent signal provides novel electron sources 
for the investigation of mesoscopic transport; 1,2 Under- 
standing their charging behavior is of crucial importance 
to assess maximum operation speed. From the theoreti- 
cal perspective such an investigation can shed light on 
the influence of quantum-mechanical processes on the 
time-dependent behavior of nanoscale systems. In this 
work, we focus on the charging of a quantum dot, which 
acts as a mesoscopic capacitor plate, the other plate is 
a macroscopic gate^^ This system, when subject to dy- 
namical modulations, is the key ingredient to realize a co- 
herent single-electron source, as shown both experimen- 
tally^ and theoretically^ The previous studies of meso- 
scopic capacitors showed that the linear low-frequency 
response regime is characterized by a quantum electro- 
chemical capacitance and a charge relaxation resistance, 
cf. the circuit in Fig. [TJ When the capacitor is at- 
tached to a single-channel lead, the charge relaxation 
resistance in the fully coherent regime is found to be 
quantized 3 -' 4 ^' 7 - demonstrating quantum mechanical non- 
locality of impedances. Theoretical discussions have ad- 
dressed the regime of weak Coulomb interaction taken 
into account on the Hartree-Fock level^ or more recently 
in the Luttinger-liquid regime, which can be treated by 
bosonization methods. ? 

In this paper, we study the time response of a single- 
level quantum dot driven out of equilibrium by either the 
fast switching (step pulse) of the gate potential or by a 
slow periodic modulation. 

The charge relaxation time after a fast switching has 
been experimentally extracted in, e.g., Ref. [l|. Also the 
time scale of spin relaxation in a single quantum dot was 
measured^ where the relaxation of the spin is due to 
spin-orbit or hyperfme interaction. For theoretical in- 
vestigations of these subjects, see, e.g., Refs. [lfj. Spin 
relaxation due to tunneling through a double dot in the 



spin blockade regime was analyzed in Ref. [Ill Here we 
are interested in the charge as well as spin relaxation in 
an interacting quantum dot, uniquely due to electrons 
tunneling to the reservoirs. 

While Coulomb interaction is weak in some of the ex- 
periments^^ side-gate coupled semiconductor quantum 
dots^ often exhibit strong local Coulomb interaction ef- 
fects and, therefore, a theory applicable to this regime 
is highly desirable. We consider a quantum dot with 
arbitrarily strong onsite Coulomb interaction, coupled 
weakly to a reservoir by a tunneling contact. We use 
a diagrammatic real-time approach to study the relax- 
ation behavior of such a mesoscopic capacitor beyond the 
Markov limit, taking into account tunneling processes 
up to next-to-leading order in the tunnel coupling (i.e. 
fourth order in the tunnel matrix elements). This ex- 
pansion permits arbitrary interaction strength U . Both 
U and temperature are assumed to be large compared 




FIG. 1: Model for the single-level quantum dot, with onsite 
Coulomb interaction U, and its circuit equivalent, (a) A step 
potential is applied to the dot at time t — to. (b) The dot 
is modulated by a slow harmonic potential, which for the 
extraction of the i?C-time we assume to be small. 
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to the line width T of the single level state. The effect 
of Coulomb interaction is manifest already in lowest or- 
der in the tunnel coupling, introducing a dependence on 
the level position of the relaxation times. We find that 
charge and spin have independent dynamics with a charge 
relaxation time that is reduced and spin relaxation time 
that is enhanced by interaction. Higher-order tunneling 
contributions lead to a difference between the RC time 
extracted from the ac admittance and the relaxation time 
extracted from the exponential decay after a fast switch. 
This is in contrast to the behavior of a classical RC cir- 
cuit and signalizes that care must be taken in extracting 
the characteristic time-scales of the system. 



II. MODEL AND FORMALISM 

We consider a quantum dot attached to a single lead 
by a tunneling contact, cf. Fig. [TJ described by the 
Hamiltonian H = H dot + H tU nnel + #lead- The single- 
particle level spacing in the dot is assumed to be larger 
than any other energy scale (temperature, Coulomb in- 
teraction, gate voltage and the highest frequency in the 
time-dependent signals) such that only one energy level 
is accessible. The dot can be described by the single-level 
Anderson model 



H. 



dot 



{t)d) a d a + Un^n x , 



(1) 



with d^dp) being the creation (annihilation) operator 
for an electron with spin a on the dot, and n a = d) a d a 
the corresponding number operator. The level position 
is assumed to depend on time via a gate potential. The 
onsite repulsion U describes the energy cost for double 
occupation and stems from Coulomb interaction (within 
the constant interaction model^). Tunneling of elec- 
trons between the dot and the lead is taken into account 
by Htunaei = J2k,a Vc l,a^a- + n - c -> where we assume a 
momentum-independent tunnel matrix element V and 
define the creation (annihilation) operators cl a (c k a ) for 
electrons with spin a and momentum k in the lead. The 
lead Hamiltonian is given by iJ lcad = J2k.a 6kC l,a c k,cr- 
We assume a constant density of states p in the leads and 
define the tunnel coupling strength r as T — 2np\V\ 2 . 
Since we are not interested in the dynamics of the non- 
interacting leads' degrees of freedom, we trace them out 
obtaining an effective description of the dot in terms of 
its reduced density matrix. We choose as basis of the 
four-dimensional Hilbert space of the reduced system the 
eigenstates, |%), of the isolated dot Hamiltonian -Hdot: 
|0) for an empty dot, \a) for a singly occupied dot with 
spin a =f, \ and \d) for a doubly occupied dot. Assum- 
ing spin-independent tunneling, we can restrict our in- 
terest to the time evolution of the diagonal elements of 
the reduced density matrix, p x . These probabilities are 
arranged in the vector P = (po->P^,Pl,Pd), whose time 



evolution is given by the generalized master equation 

dP(t) rt 



dt 



dt'W(t,t')P(t') 



(2) 



The elements W x , x > (t,t') of the matrix W(i, t') describe 
transitions from state x' &t time t' to state x a t time t. 



III. FAST SWITCHING 

We first consider the time evolution after a fast switch- 
ing event at t = to, as shown in Fig. [T] (a). The am- 
plitude of the switching is only limited insofar as we 
assume that only the single level of the dot is accessi- 
ble, i.e. the amplitude is smaller than the level spacing. 
The initial value, i.e., the equilibrium distribution before 
changing the gate voltage, is denoted by P m . In order 
to identify the relaxation time, we perform the following 
steps. First, we expand the probability vector P(t') on 
the right-hand side of Eq. © around the final time t, i.e. 

P(f) = £n=o h (f ~ tf *TeP-\t>=t. Second, we realize 
that for t' > to, the Hamiltonian is time independent, and 
hence the transition matrix elements depend only on the 
difference of the time arguments, W(t, t') — > W(t — t 1 ). 
Third, we replace the lower integration bound to by — oo. 
This is justified for times t that are larger than the char- 
acteristic time scale over which the kernel W(£) decays. 
The regime of shorter times t is not the subject of the 
present work. As a result of these three steps, we obtain 



dP{t) 
dt 



oo 1 

= E V 



w 



d n p{t) 

dt n 



(3) 



introducing the derivatives of the Laplace transforms, 



d n W := 



n 



■ f*^ dt' e-<*-*') W(t 
formal solution of Eq. ([3]) can be written as 

P(t) = cxp(Ai)P in . 



z=0+ 



The 



(4) 



The matrix A always possesses the eigenvalue 0, which 
guarantees the existence of a stationary probability dis- 
tribution. The other eigenvalues of A define the relax- 
ation rates. In the following, we concentrate on the limit 
of weak tunnel coupling between quantum dot and lead. 
This motivates an expansion of A = AW + A^ 2 ' + . . . 
in powers of the tunnel-coupling strength T, indicated 
by the superscript (i). Substituting the formal solu- 
tion for P(i), Eq. (0}, into Eq. © iteratively yields 
AC 1 ) =WW andA< 2 ) = W^+dW^-W^. Note that 
non-Mar kovian contributions, i.e. contributions propor- 
tional to <9™W, enter only starting from second order in 
r. The kernel W(z) can be computed by the real-time 
diagrammatic approach developed in Ref. Il3l . 

We are interested in the inverse relaxation times of 
charge j c = l/r c and spin 7 S = l/r s . They are found 
as the eigenvalues of A related to the left eigenvectors 
v c - (0, 1, 1, 2) - {n c ) e i (1, 1, 1, 1) and v s = (0, 1, -1, 0), 
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FIG. 2: Charge relaxation rate, 7 C , and spin relaxation rate, 
7s, in first and second order in the tunnel coupling Y for U = 
10r as a function of the dot's level position e in units of T, 
k B T = 1.5T. 



with the equilibrium occupation number (rt c ) cq . This 
yields the expressions for the time-dependent expectation 
value of the charge {n c )(t) = (0, 1, 1, 2)-P(i) and the spin 
(n s )(t) = (0,1,-1,0) • P(t) on the dot 

(n c )(t) = (n c ) oq (l - e-*/ T <=) + {n c y n e- t/T " 
(n s )(t) = (n s ) in e-^^ , 

with the initial occupation numbers indicated by the su- 
perscript "in" . The spin relaxation time can be measured 
by preparing the dot in a spin-polarized state by apply- 
ing a magnetic field and looking at the spin relaxation 
after the field has been switched off. If we retain only 
the first-order contribution in T, the inverse relaxation 
times are given by, 14 
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(5b) 



with the Fermi functions f(u>). For vanishing Coulomb 
interaction, [7 = 0, the relaxation times are independent 
of the level position e and both equal to the inverse of 
the tunnel coupling T. In the presence of a finite U, 
when the level position is in the interval —U < e < 
the dot is predominantly singly occupied. In this pa- 
rameter range, 7 C and 7 S differ strongly, even though in 
first order in the tunnel coupling a change in the spin oc- 
cupation goes always along with a change in the charge 
occupation. There are two possible states for the dot to 
relax to single occupation, namely | f), | .J,), an d con- 
sequently the charge relaxation time is decreased by a 
factor 1/2. On the other hand, in first order in tunneling 
the spin is blocked to a value different from its equilib- 
rium value due to the lack of spin-flipping processes 
leading to an enhancement of t s . This is shown in Fig. 

where the inverse relaxation times for spin, 7 S , and 
charge, j c , are plotted as a function of the new dot's level 
position after the switching, for finite U. Note that spin 



and charge have independent dynamics; this is not the 
case in the presence of a Zeeman splitting in the dot. In 
fact, in the presence of a Zeeman splitting, ^ ej,, and 
for infinite U, the first order relaxation rates are given 
by 1/n = 1 + V/(e t )/(e ; ) and l/r 2 = 1 - ^/(V)/^)- 
When we include tunneling contributions in second order 
in the tunnel coupling, we find corrections of different ori- 
gin to the relaxation times. Cotunneling, describing real 
tunneling processes via energetically forbidden interme- 
diate states, leads to coherent processes flipping the spin 
of the dot, W s f, or to transitions between zero and dou- 
ble occupation, Wod and Wdo- These terms, which agree 
with results from a standard second-order perturbation 

expansion, contain expressions (f>(e) = -^Rc'I' (j^ 

and a(e,U) = T[cf)(e + U) - 0(e)], where ^{x) is the 
digamma function and /3 = l/{k-QT) the inverse temper- 
ature, primes denote derivatives with respect to e. These 
cotunneling rates are given by 



W si = -- 

= - 



T4,"(e)+r4,"(e + U)--a'(e,U) 



2r 



e /3(2e+U) _ 1 



T<t>'{e)+T<t>'{e + U)--a{e,U) 



and Wod = exp[/3(2e+ U)} Wdo- hi addition to the cotun- 
neling terms virtual tunneling processes appear, which 
can be captured in terms of renormalization of the level 
position, e — > e + a(e,U), and the tunnel coupling, 
r^r[l + cr r (e,C/)]r^with 



<J T (e,T,U) 



T<j>'{e)+T<j>'ie + U) 

l-f(e)-f(e + U) 
l-f(e) + f(e + U) ■ 



U 



a(e,T,U) 



(G) 



These corrections are due to an interplay of tunneling 
and interaction and vanish for U — 0. We verified that 
they can also be identified in the second-order linear con- 
ductance. 

The corrections to the relaxation times read 
d_ 

'd~e 
d 



7- = ( 7(e,t/)-7| 1 )+a r ( e ,[/)7«+2W; f (7a) 



(2) 

3 

(2) 



7- = CT ( e ,[/)_7( 1 ) +(Tr ( e ,[7)7( 1 ) 

, o /(e + EQW 0d + (l-/(e))W d0 
+ l-/(e) + /(e + C0 ■ [7) 

The influence of the corrections on the relaxation times is 
shown in Fig. [2J Both relaxation rates are enhanced due 
to the higher order tunneling processes. Importantly, the 
charge and spin relaxation dynamics remain governed by 
a single time scale, each. 



IV. HARMONIC MODULATION 

We now consider a harmonic modulation of the level 
position e with frequency f2, which we assume to be much 
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smaller than T, cf. Fig. QJb). We perform a systematic 
expansion in powers of £7i 16 ' 17 The order in the expansion 
in CI, is indicated by a superscript: (i) (instantaneous) 
for zeroth order and (a) (adiabatic) for first order. While 
the instantaneous term for the average occupation of the 
dot, {n)^', immediately adjusts to the applied gate volt- 
age, its first-order correction, (n)( a ) accounts for a slight 
lagging behind. 

In the ac regime, the quantum dot can be modeled by 
its charge relaxation resistance, R q , in series with its elec- 
trochemical capacitance, We are interested in the 
frequency-dependent admittance, G(Cl), which is defined 
by the linear relation G(Cl) — I(Ct)V(Cl). Therefore we 
consider the linear regime in the gate potential control- 
ling the level position e. The expansion of the admittance 
of the dot 



G(Cl) = -iQCf, 



CI C^R q 



(8) 



for small frequencies enables to identify these quantities 
and to compute the RC time. To compare with this, we 
insert the adiabatic expansion of the occupation number 
into the charge continuity equation I(t) — —e d(n) /dt 
(Ref. [jjl) to arrive at 



1 

trc 



d(n)^/dt 



(9) 



To obtain the average dot occupation, we expand the 
master equation in powers of CI, 



= w'^P^ 



dt 



W (i )p (a ) + W^P^ + <9W (i ) 



dP« 

dt 



(10a) 
(10b) 



We refer to Ref. [Ta for the evaluation of the adiabatic 
correction of the kernel, W( a '. The instantaneous occu- 
pation probabilities, PW, and their first adiabatic cor- 
rection, p( a \ are obtained from the master equation 
Eq. (fTUf by matching orders in the tunnel coupling. We 
find that the instantaneous probabilities P^ start in ze- 
roth order in the tunnel coupling, while P( a ) starts in 
minus first order in the tunnel coupling, with p( a ' -1 ' be- 
ing proportional to the small factor Cl/T. 

A rigorous expansion in the tunnel coupling exists for 
the inverse of the -RC-time, jrc = 1/ T i?c- Close to the 
resonances, e = and e — —U, both the change in the in- 
stantaneous occupation and the adiabatic correction are 
governed by lowest order in the tunnel coupling, which 
leads to the same relaxation time as found for the charge 
relaxation time due to fast switching, 7^ = 7c 1 "*. In the 
Coulomb-blockaded regions, however, both d(ri)( l < ' /dt 
and (n)( a,_1 ^ are exponentially suppressed, and the ex- 
pression for 7^ becomes meaningless. Instead, the next- 
order corrections in T needs to be taken into account. For 
the time variation in the instantaneous occupation num- 
ber d(n)^/dt, the non-exponentially suppressed contri- 
bution to this higher-order correction is due to quantum 
fluctuations, namely cotunneling processes in which the 
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FIG. 3: Inverse RC time, 7rc, and charge relaxation rate, 7 C , 
up to second order in the tunnel coupling T for U = WT as a 
function of the dot's level position e in units of T, ksT = 1.5F 



final dot state is not changed as compared to the initial 
one. Their rates are Wo-y a ->o = Wa^o^a = — jf(j>(e) 
and W„^ a = W d ^ d = -^"(e + U). The higher- 

order correction {n)( a '°\ on the other hand, remains ex- 
ponentially suppressed. As a result, the relaxation rate 
diverges. This type of processes presumably contributes 
only to the short time dynamics after a step pulse and 
does not enter the long-time behavior described by the 
exponential relaxation. 

We find for the difference between the RC rate jrc 
and the charge relaxation rate j c 



(2) 

Trc 



(2) 
7c 



[2 - (n)^ )] <f>"{e) + (n)(*'°)0"(e + u ) 



y(l) 



d{n)^°)/de 



(11) 



Although this difference is only due to second-order rates, 
it is substantial in the Coulomb-blockade regions where 
first-order processes are exponentially suppressed. This 
is shown in Fig. [31 where the inverse i?C-time and the re- 
laxation rate after a single switching event are compared, 
taking into account tunneling processes up to second or- 
der. This difference is a consequence of the different time 
scales captured in the two expansions. Whereas in the 
response to the step we take into account the long-time 
limit but omit high frequencies, the ac response is sen- 
sitive only to the frequency at which it is tested. When 
taking [7 = this result agrees well with the respective 
limit of Ref. in the weak-coupling regime. 



V. CONCLUSION 

We calculated the relaxation times of a single-level in- 
teracting quantum dot contacted to a lead and driven 
out of equilibrium by a step-pulse potential and the RC- 
time extracted from the ac response. In the presence of 
Coulomb interaction the relaxation times depend on the 
position of the dot level, leading to different relaxation 
times for charge and spin. Corrections in second order 
of the tunnel coupling contain real cotunneling processes 



5 



and renormalization corrections to the level position and 
the coupling strength. Cotunneling processes leaving the 
dot state unchanged do not influence the exponential de- 
cay timescale leading to a fundamental difference with 
the i?C-time extracted from the admittance. 
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